Incommensurate spin-density wave in two-dimensional Hubbard model 
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Abstract. We consider the magnetic phase diagram of the two-dimensional Hubbard model on 
a square lattice. We take into account both spiral and collinear incommensurate magnetic states. 
The possibility of phase separation of spiral magnetic phases is taken into consideration as well. 
Our study shows that all the listed phases appear to be the ground state at certain parameters 
of the model. Relation of the obtained results to real materials, e.g. Cu-based high-temperature 
superconductors, is discussed. 
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Y ' Introduction 

t3 : 

Discovery of the high-temperature Cu-based superconducting compounds has given rise to in- 
vestigation of the two-dimensional (2D) single-band Hubbard model on a square lattice [T]. 
Despite its seeming simplicity this model contains the physics which underlies such complex 
electron effects as ferro- and antiferromagnetism and superconductivity. In the exactly half- 

^ ■ filling case the ground state of the model is an antiferromagnetic insulator, which agrees to 
experimental data for such mother compounds for high-T c oxides as La 2 Cu0 4 , Nd 2 Cu0 4 or 
YBa 2 Cu306- Fairly small compositional changes lead to transformation of magnetic structure 
and conducting properties of these materials. In the case of La2_ a; Sr a ;Cu04 which is the most 
extensively studied system an incommensurate magnetic order is observed by neutron experi- 
ments [2j[2]. It is generally accepted that superconducting and magnetic properties of cuprates 
are closely related. This explains a huge interest to investigation of magnetic state of the Hub- 
bard model aside from half-filling. This task appears to be a rather difficult challenge even in 

^ ' the framework of the mean-field approximation. 

Incommensurate magnetic states include spin-spiral (SS) and collinear spin-density wave 
(CSDW) states which are hardly distinguishable in experiment. In Refs. [UG3 SS structure has 
been shown to be the ground state close to half-filling and in a wide region of the magnetic 
phase diagram of the 2D Hubbard model. In Ref. [6] the possibility of phase separation has 
been studied. The main result was that the coexistence of antiferromagnetic (AF) and SS 
states is energetically favorable in the vicinity of half-filling. CSDW state is much more difficult 
for investigation than the SS one, because of need of taking into account higher harmonics 
and an accompanying charge density wave (CDW). CSDW was shown to be more energetically 
favorable than ferromagnetic (FM) and AF states [7] at certain parameters. But at the moment 
there is no data in literature concerning relative stability of SS and CSDW states. The only 
study of a so-called soliton lattice state (a set of AF domains divided by nonmagnetic domain 
walls which is similar to CSDW) shows that this state is more energetically favorable than 
the SS state at certain parameters [8]. The aim of our work is to construct a magnetic phase 
diagram of the 2D Hubbard model considering all the incommensurate magnetic states: SS 
states with the possibility of phase separation and CSDW states. 



Model and results 



The Hamiltonian of the Hubbard model: 

H = E AjJ' C j> c jV + U E n tt n U- W 

It is convenient for us to write the Hamiltonian through charge and spin densities at site: 

%% = ^j-f' (2) 

We make the mean-field approximation which consists in the neglect of the terms of second 
order of smallness in the expansion of the charge and spin densities in a mean value and a 
fluctuation: 

nj^-nj + 2njnj (3) 

Sf -Sf + 2S?S?. (4) 

For consideration of CSDW we take a unit cell which size L is multiple of the CSDW period. 
If the wave vector is directed along the diagonal of lattice Q = (Q, Q) than a unit cell is chosen 
elongated in that direction as well. Taking into account the periodicity of magnetic structure, 
we have S? t = Sf, n- iy i = ni, where I is the site number in the unit cell. 

Then the Hamiltonian takes the form: 

H = K + V, (5) 
where the kinetic energy (only nearest neighbor hopping is considered): 

^=*EEE c jW>- ( 6 ) 
jj' i,v ° 

and the interaction: 

V = UN Y. (Sf - % ) + WY. {%,, - S?S|,) . (7) 

The Fourier transformations lead to the following form of the Hamiltonian: 

K = 1 E + + <l,*CKi-iAe~ ik * + e~ lky )) (8) 

V = UNj2(sl 2 -^)+Uj2[(j- S?) 4 M c kM + (I + Sf) < u c k , u ] . (9) 

l ^ ' k,l 

Similarly one can obtain the Hamiltonian for parallel states (Q,7r) and (0, Q). There con- 
sidered are the diagonal and parallel CSDW states only because of their relative simplicity 
and on the basis of data obtained for SS states for which only diagonal and parallel states 
were shown to be the ground ones [5]. To diagonalize this Hamiltonian we should diagonalize 
2L x 2L matrices. The total energy of the system can be find by solving the self-consistent 
system of equations on charge and spin densities at sites. We suppose that these densities can 
take any values through the self-consistency procedure, i.e. we do not fix a profile of CSDW 
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Fig. 1: (Color online) Magnetic phase diagram of the 2D Hubbard model on a square lattice. 
Bold green lines denote calculated boundaries of CSDW state, blue lines denote phase transi- 
tions between homogeneous SS states, red lines denote phase transitions between separated SS 
states. Regions of CSDW state are dark shaded (green), regions of phase separation are light 
shaded (grey). 




Fig. 2: The spatial profiles of spin (lower plot) and charge (upper plot) densities as a function 
of lattice sites for U = 8, n = 0.8. 

and CDW. For fixed U / 1 and n — \ Yl\ ^i we §° over different wave vectors Q and search 
one that provides a minimum of the total energy. We compare this minimal value with the 
energy of the ground SS state with phase separations which is calculated by the numerical 
scheme described in Ref. [6|. Thus we find the overall ground state considering different in- 
commensurate magnetic phases. The results of calculations are presented in the phase diagram 
(Fig. [p. We see that the CSDW states appear to be more stable than the SS ones for small 



and intermediate U/t both in the vicinity of half-filling and far from it. It should be noted 
that the wave vector which provides the minimum of the energy can be, as a rule, calculated 
by formula Q = nn (for (Q,Q) and (Q,tv) states). The boundary of CSDW and PM states 
is supposed to coincide with the boundary of SS and PM states — accuracy of calculations 
does not allow to determine this line more precisely. On Fig. [2] spatial profiles of CSDW and 
CDW are presented for U/t = 8, n — 0.8, Q = (0.87T, 0.87r). One can see a strong deviation 
from the harmonic wave form. The larger U / 1 the closer magnetic structure to a soliton lattice 
state, i.e. antiferromagnetic domains separated by nonmagnetic domain walls. As follows from 
the figure, CDW has sufficiently large amplitude — from 0.47 at nonmagnetic sites to 0.96 on 
sites with maximal magnetization. Such a strong charge inhomogeneity looks rather unrealistic. 
An additional Coulomb energy caused by this inhomogeneity is not taken into account in the 
framework of our model, but it is rather obvious that it should substantially raise the energy 
of the CSDW so that it would not be the ground state. 

Summary 

We have constructed the magnetic phase diagram of the 2D Hubbard model in the mean- 
field approximation. The calculations show that the CSDW state is the ground state in a 
wide region of the model parameters. It is more stable than the uniform SS state and the SS 
state taking into account phase separations. We did not take into consideration an additional 
Coulomb interaction rising from the charge inhomogeneity of the system. We believe that this 
interaction could be sufficiently strong to suppress CSDW state. That is why we think that the 
true nature of the CSDW state stabilizing in real systems is beyond the mean-field treatment 
and is still open to question. 
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